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We investigate the structure of magnetic domain walls in a classical double exchange ferromagnet, evaluating 
domain wall energies and charges. Three different cases are studied: (i) a conventional smooth Bloch wall, 
(ii) an abrupt Ising-type wall, which is shown to have lower energy at small values of carrier concentration, 
and (iii) stripe wall, corresponding to the two ferromagnetic domains being separated by a stripe of another, 
antiferromagnetic, phase. General aspects of energy balance and geometry of phase-separated states are 
discussed in this context. It is speculated that domain walls of the latter type may be responsible for the 
unusual transport properties of certain manganate films. 
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I. INTRODUCTION 

The unusual micro-magnetic properties of colossal mag- 
netoresistance (CMR) compounds are presently subject 
to intensive experimental investigation [1-11]. In these 
studies, special attention is paid to the interplay between 
magnetic domain structure and transport properties of 
the system. Aside from possible technological applica- 
tions (associated with the large low-field magnetoresis- 
tance [6]), the strong effect of magnetic domain walls 
on conduction properties, as found in strained epitaxial 
films of La .7Ca . 3 MnO3 (Refs. [5,6]), Pr 2/3 Sr 1/3 Mn03 
(Ref. [6]), and Lao.7Sr .3Mn03 (Refs. [6,7]), raises a gen- 
uine physical problem. Indeed, given the relatively small 
expected value of the easy-axis magnetic anisotropy, 
the usual Bloch (or Neel) domain wall would be rather 
smooth and broad. Thus, carrier scattering off the Bloch 
walls could not appreciably affect transport properties of 
the system. The measurement of magnetic domain walls 
contribution to the resistivity therefore leads to the con- 
clusion [12] that the domain walls arising in the samples 
studied in Refs. [5-7] have an unusual, non-Bloch struc- 
ture. It has even been suggested [6] that the double ex- 
change interaction, which is responsible for the ferromag- 
netism of doped manganese oxides, cannot possibly ac- 
count for such poorly-conducting magnetic domain walls. 
While the origins of this suggestion may be traced to 
the widespread but ill-founded notion that the magnetic 
properties of double exchange systems can be adequately 
described by an effective Heisenberg model, the peculiar 
physics of domain walls in double exchange ferromagnets 
has not yet been addressed theoretically. 



In the present article, we consider the standard single- 
orbital double exchange model with the following Hamil- 
tonian: 
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Here Cj a (with a =|, J) are the electron annihilation op- 
erators, and the vector a a/3 is composed of Pauli ma- 
trices. Jh is the strength of Hund's rule ferromagnetic 
coupling between the spins of carriers and the core spins 
Si, which also interact with each other via the direct anti- 
ferromagnetic Heisenberg exchange J. The core spins are 
assumed to be classical (S ^> 1), and the easy-axis single- 
ion anisotropy Kj S 2 is included in order to account for 
the finite Bloch wall energy. The lattice is assumed to 
be square, which is thought to be more appropriate than 
the three-dimensional cubic one for modelling the thin 
films studied experimentally; the extension of our anal- 
ysis to the three-dimensional case is straightforward but 
cumbersome, and is expected to yield similar conclusions. 
The electron spectrum in the ferromagnetic state is given 
by [13] el' 1 = e^J H /2 with eg = -t(cos k x +cos k 2 ). We 
consider the experimentally relevant half-metallic case, 
when owing to a sufficiently large value of Jh, the car- 
rier band in the ferromagnetic phase is completely spin- 
polarised. Thus the value of chemical potential, denoted 
(i — Jh /2, must lie below the bottom of the spin-up sub- 
band, /i < Jh — 2. We note that t in Eq. (1) corresponds 
to 2t in a different notation sometimes used elsewhere in 
the literature; it should also be pointed out that below, 
the conduction electron (rather than hole) density is de- 
noted by x. Throughout the paper we use units in which 
hopping t and the lattice spacing are equal to unity, and 



1 



we consider the zero-temperature (T = 0) case. 

Below we consider domain walls of three different 
types, which are relevant for different values of pa- 
rameters characterizing the double exchange magnet at 
low temperatures. These are conventional Bloch walls, 
abrupt (Ising) walls, and stripe walls, formed by a stripe 
of antiferromagnetic phase inserted between the two fer- 
romagnetic domains. 

We begin in Sect. II with the usual smooth Bloch wall. 
The Bloch wall energy depends on spin stiffness D and 
anisotropy strength in a usual way [14], reflecting the 
fact that the long-wavelength properties of double ex- 
change ferromagnets are adequately captured within an 
effective Heisenberg description (cf. Ref. [15]). In double 
exchange systems, Bloch walls carry an electric charge, 
which we also evaluate. Our results suggest that mag- 
netic domain walls arising in homogeneous (single-phase) 
double exchange ferromagnets at the intermediate dop- 
ing levels typically have Bloch structure, and therefore 
cannot significantly affect the resistance of the sample. 

Within the context of recent experiments, the possi- 
bility of sharp changes in magnetisation direction within 
a domain wall is of particular relevance. This scenario 
has been discussed for a long time (see, e.g., Ref. [16]) 
and it is important to consider it in some detail. There- 
fore in Sect. Ill we treat the extreme case of an abrupt 
(Ising-type) domain wall, when the sign of magnetisation 
is reversed over one lattice link. A wall of this type, which 
in the T = 0, S — > oo limit is impenetrable for carriers, 
would strongly affect the transport properties of the sys- 
tem. The energy cost of an abrupt wall originates from 
the underlying non-perturbative scattering problem for 
conduction electrons. The corresponding physics is thus 
completely non-Hcisenberg. We derive expressions for 
energies and charges of abrupt domain walls running in 
two different directions (along a crystal axis and diag- 
onally), and for all values of the Hund's rule exchange 
constant, Jr- While for small values of carrier density, 
x <C 1, the energy of an abrupt wall is lower than that 
of a Bloch wall (which may be relevant for certain mag- 
netic semiconductors), this does not generally hold at 
the intermediate doping levels. In the latter case, an 
abrupt domain wall is preferred only for very large val- 
ues of anisotropy K ~ DS, or for the case of very finely 
tuned parameter values, providing for an almost exact 
balance between the ferro- and antiferromagnetic tenden- 
cies of the system. It would be unrealistic to expect that 
such a fine-tuning (within one per cent in the values of 
J, Jh, and x in a single-phase system) can be achieved 
by different experimental groups in a reproducible way. 

In addition, it also turns out that these parameter val- 
ues typically correspond to the system being unstable 
with respect to phase separation. As explained in Sect. 
IV, the latter phenomenon has a double effect: (i) the 
carrier density within the bulk of the ferromagnet is now 
determined by the condition that the thermodynamic po- 



tentials of the two phases must be equal to each other; 
this condition effectively pins the parameters of double 
exchange ferromagnet in the region where the energy of 
an abrupt domain wall is relatively low. (ii) Energy of 
an abrupt domain wall can be further lowered by insert- 
ing a stripe of antiferromagnetic phase between the two 
ferromagnetic domains. Since the two phases are charac- 
terized by different values of charge density, one cannot 
treat this situation properly without taking into account 
the effects of Coulomb interaction. We use a somewhat 
simplified treatment to estimate the energy and width of 
a stripe domain wall. It turns out that within a certain 
range of parameter values, the energy of a stripe wall 
can be lower than that of a Bloch wall, so that magnetic 
domain walls in a phase-separated system are actually 
of the stripe type. In particular, this situation is real- 
ized when the antiferromagnetic phase occupies an ap- 
preciable area of the sample (of the order of 15 % of the 
net area, or possibly more), provided that the easy-axis 
anisotropy constant K is not too small. Due to insu- 
lating properties of the antiferromagnetic phase, carrier 
transport across the stripe wall is strongly suppressed, 
leading to a substantial domain wall contribution to the 
sample resistance. On the other hand, ferromagnetic area 
within a single magnetic domain remains well-connected, 
and phase separation is therefore not expected to sig- 
nificantly affect the intra- domain metallic conductivity. 
Analysis of the data of Refs. [5,6] reveals a correlation 
between the film thickness, dielectric properties of the 
substrate, and the appearance of domain wall resistance, 
which seems to agree with anticipated conditions for the 
stabilisation of the stripe walls. 

Details of calculations are relegated to the Appen- 
dices, which also include a brief discussion of the three- 
dimensional case. 

The relevance of our findings in the context of recent 
experiments on manganate films is further discussed in 
Sect. V. We suggest that the domain walls observed indi- 
rectly in the transport measurements of Refs. [5-7], and 
directly in Ref. [8] are in fact the stripe walls, introduced 
in Sect. IV. 

II. BLOCH WALL 

The structure of domain walls in conventional Heisen- 
berg ferromagnets has been understood long ago [14]. 
These are smooth, long- wavelength Bloch walls [17], and 
their surface tension (energy per unit length) Sb and 
width Ib are determined by the spin stiffness D of the 
system: 

S B = ■ (DS) , l B = ^DS/K . (2) 

Since the unusual transport properties of the domain 
walls are found only in certain strained films at a spe- 
cific doping level [5-7], we expect that in most cases, 
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domain walls in the CMR materials also have Bloch-like 
structure. We will now study the relationship between 
the properties of Bloch walls and the parameters of our 
model Hamiltonian, Eq. (1). 

The appropriate value of D can be extracted from the 
known spin wave spectrum of a classical double exchange 
ferromagnet [18] (see also Ref. [15]): 
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Here, A*" is the number of lattice sites, and is the 
Fermi distribution function. Eq. (3) is valid for any 
dimensionality d and for an arbitrary electron dispersion 
law (with el'* = T-Jff/2 + eg). For the case of the 2D 
tight-binding model (1), we obtain: 
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Here and below, Yj and Y" 2 denote the following complete 
elliptic integrals: 



(5) 



= de^/dk is the electron velocity, and the kinetic en- 
ergy of the band is given by 



(6) 



Note that because of the numerical prefactor entering 
Eq. (4) , the value of D is at least an order of magnitude 
smaller than that of the band energy, E. 

At low doping level, x -C 1 ~ Jh, Eq. (4) yields 



DS = -J+- 4 x--nx -—, 

whereas at half- filling, x = 1, we obtain 
DS = —J — 1/(4 Jh) . 



(7) 
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The second term in Eq. (7), which is proportional to 
the band energy (E w — 2a; at low x) represents the 
leading-order double-exchange (ferromagnetic) contribu- 
tion. The last terms in Eqs. (7-8) indicate that the 
effect of finite Jh (as opposed to Jh — » oo) is similar to 
that of an increase in the value of direct superexchange, 



J. This conclusion is justified physically, since at finite 
Jh an effective antiferromagnetic interaction arises due 
to virtual transitions between the two components of the 
spin-split band much like a usual superexchange, which 
is due to transitions between different bands. Below we 
will see how this qualitative analogy [19] manifests itself 
in other properties of the system - its validity is clearly 
not restricted to the spin stiffness evaluation. This in 
turn suggests that many of the features of (more compli- 
cated) finite- Jh systems can be modelled by treating the 
Jh — * oo case with an appropriately increased J. 

The doping dependence of spin stiffness for three differ- 
ent values of Jh (Jh — * oo, Jh — 8, and Jh = 4 for solid, 
dashed, and dashed-dotted lines, respectively) and J = 
is shown in Fig. 1 (a). For the case of finite Jh, the com- 
petition between effective antifcrromagnctism and double 
exchange-induced ferromagnetism, taking place at suf- 
ficiently small 1 — x, is resolved via phase separation 
[20-23]. This means that the homogeneous ferromag- 
netic state becomes thermodynamically unstable as the 
electron concentration x exceeds certain critical value. In 
Fig. 1 (a), the values of DS within the respective ther- 
modynamically unstable regions are plotted with dotted 
lines. When the superexchange J > is present, this crit- 
ical value, which depends also on Jh, decreases further. 
In addition, another region of phase-separation instabil- 
ities arises at low electron densities [21,22]. 

Within a Bloch wall, misalignment of the neighbouring 
ionic spins leads to a renormalisation of carrier hopping 
coefficient [24]. Indeed, the Hamiltonian (1) can be re- 
written in terms of new fermions (and whose 
spin is aligned (antialigned) with the classical ionic spin 
Si at the same site: 



n = -\ E (tf4 a d 3 , + tfd] a d t0 



(9) 



+ Jj t E K*i - 4 t *0 + i E *4 - J* E (^) 2 • 



(i,3) 



Here, the matrix t a/3 is given by 



a/3 = / CiCj + jW-MSiSj 

\ -e^'SiCj+e^CiSj 
9, 



-e-^dSj+e-^SiCj 
CiCj +e i (*«-*')S i S j 



are the polar co-ordinates of the spin Si 



Ci = cos , S t = sin -± , (10) 
and 6 

In the bulk of the ferromagnetic state, reduces to 
a unit matrix, but inside the domain walls, the values 
of both diagonal and off-diagonal elements are changed. 
Thus, the bandstructure (and hence the carrier density) 
within the wall differs from that in the bulk, and we come 
to the conclusion that Bloch walls are charged. We will 
now evaluate the surface charge (Tb of a Bloch wall in a 
double exchange ferromagnet. 
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Let us suppose that the Bloch wall runs along the [11] 
direction of the lattice diagonal, and choose the y axis 
to be perpendicular to the wall. We also choose the co- 
ordinates in spin-space in such a way that <pi = 0, and 
note that 6i does not depend on x. In other words, the 
spin configuration is composed of ferromagnetically or- 
dered chains running in the x direction, with the inter- 
chain and intrachain distances given by l/y/2 and a/2 
respectively. It is then convenient to Fourier-transform 
the fcrmion operators in the x direction only, according 
to 

/ 2 \ 1/4 

d a (x,y)=l-) Y,e ik * x ^d a (k x ,y) , \k x \ < tt . 



(11) 



Then the first two terms in Eq. (9) can be re-written in 
the form 

K=-J2 C0S T t a ^y,y + ^)dl(k x ,y)d fj (k x ,y + 



+h.c. 



y,k x ,a,/3 

J 
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+ -^- [d\(k x ,y)d l (k x ,y) - d\(k x ,y)d^{k x ,y) 



(12) 



n = x 

r-l 



5x 



(16) 



2de x 



n/2 n^/l-el 7-2 



{Py(e y ) + 5p y (e x ,e y )}de y . 



Here, p y = l/(7r^/4 — e y ) is the value of density of states 
at fixed k x in the ferromagnetic state, and 8p v is the 
correction arising at dO/dy ^ 0. Then the change in the 
carrier density due to a non-zero value of (89/ dy) <C 1 
can be evaluated (to leading order in Se y ) as 



Sx 



<K(~) 



1 2de x 

H/2 TTy/T— e 

v/e x 



Sn y (—) 



6p v (e x , € y )de y w -Se y p y ( — ) . (17) 



Using Eq. (15), we obtain after some algebra 
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Finally, given the known profile of 9(y) in a Bloch wall 
[14], cos9(y) — tanh(y/Zs), we find the following expres- 
sion for the charge of a Bloch wall per unit length: 



which we will also use in Sect. Ill below. 

In the ferromagnetic state, the subsequent Fourier 
transformation in the y direction according to 
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d a (k x ,y)=^—j £e ifc ^ Q , \k y \ 



<ir. (13) 



yields the spectrum, 
Jh 



4 4 



-2 cos — cos k y . (14) 



The variation of spin direction within a Bloch wall cor- 
responds to the long- wavelength limit, Ib 3> 1, of con- 
tinuum micromagnetic theory. Then one can define 9(y) 
as a continuous function, and the angle formed by the 
spins S(x,y) and S(x',y + 1/V2) on the neighbouring 
chains is given by (d9/dy)/V2. For the case of a con- 
stant value of 86 jdy <C 1, the spin-up fermion spectrum, 
et = — (Jjj/2) is obtained from Eqs. (10-12) (upon 
Fourier transformation, Eq. (13)). When 89 jdy ^ 0, the 
quantity eg is only approximately factorisablc, 



eg = cos -y- (e y (k y ) + Se y (k x , k y )) , e y = -2cosfc y , 
4 k 

5e y = (89/8y) 2 (cos k y — cos sin 2 k y )/8 . (15) 

Jh 2 

The value of carrier density at a fixed value of chemical 
potential is then given by [25] 



where e is the absolute value of electron charge. In eval- 
uating c b as — e J Sxdy, we used the adiabatic approx- 
imation, which is valid in the long-wavelength limit of 
Ib S> 1- As expected, a similar calculation for a Bloch 
wall running parallel to a lattice direction yields the same 
result (19): Bloch walls have a well-defined continuum 
limit, and both their energy [26] and charge are indepen- 
dent of the orientation on a square lattice. 

We note that at J > 0, the Bloch wall can be stable 
only as long as the chemical potential at the centre of 
the wall (where the band-narrowing effect is most pro- 
nounced) lies above the bottom of the carrier band. In 
other words, the value of x+5x with Sx given by Eq. (18) 
should remain positive at y = (otherwise, there would 
be no carriers and hence no carrier-mediated ferromag- 
netic interaction near the centre of the wall) [27]. Since 
the wall is smooth, Ib ~> 1, this condition is important 
only at the low-doping limit of x <C 1, when it reads [28] 



l%i:xDS > K. 



(20) 



This is clearly violated at sufficiently low x. We will see 
that in this case the domain wall is in fact abrupt [Sect. 
Ill, Eq.(21)]. 

According to Eq. (19), the charge of the Bloch wall, 
which is inversely proportional to its width, decreases 
with decreasing anisotropy strength: as oc \[K. At 
small values of electron density x <C 1, Jh, wc find 
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(7b = e/(87rZs). The behaviour of <jb at the intermediate 
doping levels can be inferred from Fig. 1 (b), where the 
quantity C(x) (see Eqs. (18-19)) is plotted for different 
values of Jh- We suggest that the experimental deter- 
mination of <jb may help to distinguish Bloch walls from 
abrupt or stripe domain walls (see Sections III-IV be- 
low), which typically carry larger charge. On the theory 
side, the effect of Bloch wall charge on the carrier trans- 
port across the wall should be considered. 

Throughout this section, we assumed [29] that the 
Debye-Hiickel screening radius is large in comparison to 
Ib- This appears to be plausible, especially in view of rel- 
atively large values of dielectric constants, characteristic 
for the highly-polarisable oxides. We will briefly discuss 
the magnitude of Coulomb correction to the Bloch wall 
energy, Sb, in Appendix B [Eq. (B9)]. In the opposite 
case of strong screening, the charge of a Bloch wall will 
vanish. 

III. ABRUPT WALL 

The appreciable contribution of magnetic domain 
walls to resistivity, as observed in certain ferromagnetic 
strained CMR films [5-7] , suggests the possibility of non- 
Bloch walls arising in these systems. Indeed, in order to 
scatter the carriers effectively domain wall must have a 
non-smooth structure, characterized by abrupt changes 
in spin direction. An abrupt (Ising-type) domain wall, 
shown in Fig. 2, represents an extreme example of such 
a structure. 

Unlike the Bloch wall, abrupt wall represents a lattice 
problem (as opposed to a long- wavelength one). There- 
fore the properties of an abrupt wall depend on its orien- 
tation with respect to the lattice, and one has to distin- 
guish between, e.g., diagonal (Fig. 2 a) and vertical (Fig. 
2 b) walls. We note that a similar feature would also arise 
for domain walls in an Ising fcrromagnct - indeed, the 
number of cut ferromagnetic links per unit wall length 
is different for vertical and diagonal walls. In a classi- 
cal double exchange ferromagnct, the standard double 
exchange mechanism forbids carrier hopping across the 
abrupt domain wall [30]. Owing to the anisotropy of 
the carrier spectrum (as manifested in a non-spherical 
shape of the Fermi surface), the carrier contribution to 
the abrupt wall energy is again orientation-dependent. 

In order to show that abrupt domain walls can actually 
arise in double exchange ferromagncts, we will first turn 
to the low-doping limit, assuming also that Jh = 

oo and J = 0. Since the Fermi momentum is small, 
Pp = A-kx <C 1, carrier dispersion can be approximated 
by the free-particle dispersion law, el w const + (k 2 /2). 
The energy of an abrupt wall is therefore equal to that 
of a partition inserted into an ideal spin-polarised Fermi 
gas, which can be easily estimated. 



Let the ideal Fermi gas be contained in a rectangular 
box of the size L x x L y . According to the uncertainty 
principle (or alternatively to the usual rules of momen- 
tum quantisation), the difference between the allowed 
values p y of the y-component of momentum can be esti- 
mated as Sp y ~ 1/L y . Suppose now that a flat partition 
perpendicular to the y axis has been introduced, divid- 
ing the box in half. This shifts each allowed momentum 
value: p y ^ — > p$ + 5p y ^ with \5py^\ <~ 5p y . The signs of 
Sp^ are chosen in such a way that the energy shift of each 
individual electron level is positive: Se(p x ,p y ) ~ \p y \5p y . 
The net energy change associated with the partition is 
thus given by L x L y J np\p y \d 2 p/ L y , or <~ x 3 / 2 per unit 
length of partition [31]. 

Thus, we find that the energy of abrupt domain wall 
in a double exchange fcrromagnct is given by Sa ~ x 3 ^ 2 . 
The numerical coefficient can be obtained by an exact 
treatment [see below and Appendix A, Eqs. (A16-A17)], 
yielding Sa ~ 4y / 7ra; 3 / 2 /3. Comparing this with the 
Bloch wall energy, Sb ~ \[Kx (see Eq.(2)), we find that 
the abrupt wall energy is lower, Sa < Sb, as long as 

x 2 < 9K/16ir. (21) 

We note that according to Eq. (20), Bloch walls become 
altogether unstable at x 2 < K/An. 

It appears to be very difficult to rigorously address the 
question whether in the region specified by inequality 
(21) the abrupt wall actually represents the optimal spin 
configuration. We are, however, able to verify [see Ap- 
pendix A, Eqs. (A16-A17)] that as long as x 2 < K/ir, 
the abrupt domain wall is stable with respect to small 
"smearing" perturbations (shown schematically in Fig. 
3) involving spins adjacent to the domain wall on both 
sides. This provides a strong, albeit variational, argu- 
ment for the overall stability of abrupt walls. 

We now turn to exact calculation of energies and 
charges of abrupt walls for all values of x, Jh, and J, 
beginning with the evaluation of the electronic contribu- 
tion to the energy of an abrupt diagonal wall. 

Following the Fourier transform, Eq. (11), the elec- 
tronic terms in the Hamiltonian of the uniform ferro- 
magnetic phase take the form (cf. Eq. (12)) 

n kx = -§ ]r{4(fc.,y + ^=H(fc^)+ 

y I V 

+ d\{k x ,y+ -j=)di(h x ,y) +h.c.| + 

+^^[d\(k x ,y)d l (k x ,y) - d\(k x ,y)d ] (k x ,y)^ , (23) 
y 

where Q — 2 cos(fc x /2). The abrupt diagonal domain 
wall parallel to the x axis results in a perturbation of the 
Hamiltonian (23), Hk^ — » Hk x + V klc , with 
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+ 



x (1 - cos V>) + |dJ T c?i| + c^dix | (1 - sin 2^>) 

+ |rfJ T c?ij. - ^o^it} cos 2^ + |d^ 1T d x - d^_ 1 ^dof + 

+d\^d 2 i - 4i rf 2t} sin^ + h.c. . (24) 

Here we denoted d a (k x ,i/y/2) by d iQ and allowed for a 
smearing perturbation, ip <C 1, as shown in Fig. 3. It 
is convenient to re-write the operator Vk x in a diagonal 
form, 



Vk x = ^ Aia\ai , a\aj + Oj-aJ = <5 



(25) 



i=i 



Expressions for both the eigenvalues Ai and the operators 
at are given in Appendix A. 

In the absence of a domain wall, the electronic contri- 
bution to thermodynamic potential of a double exchange 
ferromagnet at a temperature T can be evaluated as 



= 



2irV2 J 2^ rv < 



/L x dk x f 
2ttV2 J 



dev tot (e,Q)if(e)de, 



2ttV2 

= — Q cos k y , (/3(e) = —Tin 



exp 



(26) 



Here, L T and L,, are the dimensions of the sam- 



ple, vtot{t,Q) = L y ^/~2/(ny/Q 2 — e 2 ) is the total den- 
sity of states at a fixed value of Q [i.e., with k x — 
±2arccos(Q/2)], and the factors \[2 originate in mo- 
menta rescaling implied in Eqs. (11) and (13). 

When the domain wall perpendicular to the y axis is 
introduced, the associated perturbation Vfe x , Eqs. (24- 
25), gives rise to a correction [32] in the density of 
states, u t ot(e, Q) — > v to t{e, Q) + Si>(e, Q). Introducing the 
Lifshits-Krein spectral shift function [33] £(e, Q) accord- 
ing to 8v = —d£/de, we find for the electronic contribu- 
tion to the domain wall energy, 

^ de5v(e, Q)ip(e) = 



2irV2 
dh 



r dk x r 

J 2irV2~J 



de*(e,Q)/(e) 



(27) 



Here, the zero-temperature value for the Fermi distribu- 
tion function, /(e) = 9(p — e), can be substituted. 

For a given value of k x , the operator T4 x represents 
a local perturbation of a one-dimensional Hamiltonian 
7ife x . Thus, the dependence of £ on Q is only parametric 
[32] , and the value of £ can be found from the standard 
formula [33] (see also Ref. [34]): 



£(e,Q(fc x )) ArgDet 



l-G(e--J H -iO,Q)V k . 



(28) 



where G(£, Q) = (£• 1 — Hk^ 1 is the resolvent operator 
at a given value of k x , and 1 is the identity operator. In 
the basis containing the states a||0) ( where |0) is the 
vacuum state), the determinant on the r. h. s. of Eq. 
(28) is that of an 8 x 8 matrix, 8ij — MijAj, with 



Mij = E / 



dfc B V2(0|o«|fcg)<fcg|at| 



<*=t,r 



2nQ E a + cos k y — iO 



(29) 



Here 



E x =e/Q, E l = (e - J H )/Q , (30) 
and \ky) are properly normalised Bloch wave states, 

l fc y> = ^E c ~ ife ^<( fc -y)l°>' ( 31 ) 



(k«\k>/) = 2nS(k y ~k> y )5 af} . 



After a straightforward, if somewhat laborious, calcu- 
lation we obtain 



£(e,Q)=£ (0) (e,Q) + <5£, tg7r£ 



(o) 



E-[Ei - 1 



A 72 

ttQ z 



l-E. 



E 2 - 1 



(32) 



The final expression for the energy of an abrupt diagonal 
domain wall per unit length is then given by the trace 
formula, Eq. (27), with additional contributions from 
direct superexchange and single-ion anisotropy: 

S d = S d 0) + Z d 4, 2 = -2V2J + 2V2(2 J + K)yj 2 + 



QdQ 



■■ ^(QEt,Q)8(n - QEf) . (33) 



The energy of a vertical abrupt domain wall is calculated 
very similarly (see Appendix A), yielding the result 



S v = Sf> + Z v ^ 2 = -2 J + 4(J + K)i> 2 + 
+ l_f_de, 
k J-i 



de 2 i{e 2 )6(fj, -e\- e 2 ) . 



(34) 



Here, £(e 2 ) is equal to £(e 2 , Q) as given by Eqs. (32) with 
Ei = e 2 , Ei = e 2 - Jh, and Q=l. 

The spectral shift function, Eq. (32), also contains 
information about the abrupt domain wall charges. In- 
deed, spectral shift function £(e) generally measures the 
number of energy levels that cross the given energy value 
e as a result of a perturbation. Thus, the change in elec- 
tron density at a fixed value of k x is given by — £(//, Q), 
yielding the charge of an unperturbed (tp — 0) abrupt 
diagonal wall: 



G 



V2 . , 
Od = e / £ 



(o) a 



dQ 



For a vertical wall, we likewise obtain 



<>■„ = -/ £ (0) - eisgn/x) 4^=? ■ 



(35) 



(36) 



Here the function £(e) is defined in the same way as in 
Eq. (34) above. 

We have conducted a thorough numerical investiga- 
tion of Eqs. (33-34). Doping dependence of the abrupt 
wall energies for different values of Jh is illustrated in 
Fig. 4 (a). Comparing these with Fig. 1 (a), we con- 
clude that at the intermediate doping levels, abrupt wall 
energies are typically several times larger than spin stiff- 
ness, DS. Therefore in the physically relevant case of 
small anisotropies, K <C DS, Bloch walls will typically 
have a significantly lower energy [see Eq. (2)]. We note 
that including antifcrromagnetic superexchange, J > 0, 
would lead to a decrease in S d °^ relative to si°^ [as fol- 
lows form Eqs. (33-34)]. In particular, this can yield [35] 
S d °^ < si°^ at small values of x. 

The charges, ad and a v , of the abrupt domain walls are 
plotted in Fig. 4 (b). We see that at the intermediate 
doping values, the electric charge per unit length is of the 
order of O.le, in a marked difference from weakly-charged 
Bloch walls [cf. Fig. 1 and Eq. (19)]. 

With increasing antifcrromagnetic interactions in the 
system (that is, either with increasing J or with decreas- 
ing Jh) spin stiffness, as well as the abrupt wall energies, 
will eventually change sign. Near this point, there might 
be a region where DS is still positive, while either si°^ 
or sjp is smaller than the Bloch wall energy, Sb- This 
is due to the fact that Sb and S^ d depend on J and 
J H in different ways. Such a situation is illustrated in 
Fig. 4 (c), showing the domain wall energies in a dou- 
ble exchange ferromagnet with x = 0.55 and Jh = 4 as 
functions of J. The solid line corresponds to the Bloch 
wall energy Sb, whereas the the vertical abrupt wall en- 
ergy S^ , is represented by a dashed line. The value 
of easy-axis anisotropy constant, K, is varied with J in 
such a way that Bloch wall width, Is [see Eq. (2)], is 
always equal to 5. We see that Sb > si ' 1 > for 
0.0143 < J < 0.0148. Since the quantity Z v [see Eq. 
(34)], represented by the dotted line is positive through- 
out the si°^ < Sb region, one is tempted to conclude 
that the abrupt wall is indeed stable in this region. How- 
ever, it is easy to verify that for J > 0.0107, the homoge- 
neous ferromagnetic ground state becomes unstable with 
respect to phase separation into ferro- and antiferromag- 
netic regions. It appears that this represents the general 
situation, i.e., that at the intermediate doping range the 
inequality Sb > S^ (or Sb > S d ^) cannot be satisfied 
within the thermodynamically stable region. In Sect. IV 



below, we will argue that the phenomenon of phase sep- 
aration can affect the magnetic domain wall structure in 
a profound way. Here we merely note that even if phase 
separation is suppressed due to some mechanism (e. g., 
enforcing electric neutrality on the microscopic level), the 
parameter region where either S v or Sd is smaller than 
Sb (but the stiffness D is still positive) would still be 
very narrow, requiring one to fine-tune the values of J, 
Jh, K, and x to within a fraction of a per cent [36]. It 
is therefore very unlikely that such a situation can be 
realized experimentally in a reproducible way. 

Expressions (33-36) can be further simplified in the 
limiting cases of small carrier density, x <C 1, or large 
Hund's rule coupling, Jh — > oo (see Appendix A). Ex- 
pressions (A12-A15), valid in the Jh — > oo limit, can be 
used to estimate the values of Sd and S v at sufficiently 
large Jh throughout the entire range of dopant concen- 
trations. 

As discussed in the beginning of this section [see Eq. 
(21)] , the domain walls become abrupt at the low-doping 
limit of x <C 1. In this case, the abrupt wall energies and 
charges are given by Eqs. (A16-A19). The doping depen- 
dence of domain wall energies in this region is illustrated 
in Fig. 5. The value of Hund's rule coupling is taken to 
be Jh — 0.1, and K(x) = D(x)S/25 again ensuring that 
Ib = 5. We see that the the energy of an abrupt vertical 
wall (dashed line) is lower than that of a Bloch wall (solid 
line), S ( v 0) < S B , for all x < 0.0027, and the stability of 
abrupt domain wall is further evidenced by the fact that 
the the quantity Z v (dashed-dotted line) is positive for 
x < 0.0063. For this choice of parameters [37], the value 
of S^ is just above that of Sv°\ and we find S d °^ < Sb 
&tx< 0.0026, and Z d > for all x < 0.0074. Bloch wall 
becomes unstable [see Eq. (20)] at x < 0.0008 (dotted 
line). We note that allowing for a larger value of K would 
have broadened the region where abrupt walls have lower 
energy; however, Eqs. (2) are valid only in the !b > 1 
case. 

The data shown in Fig. 5 are for a system with no 
direct superexchange, J = 0; including a small J > 
would give rise to a phase-separation instability at small 
x [22], which may or may not cover the entire region 
of S^,°] < Sb- While no study of domain structure in 
the electron-doped manganates has been reported so far, 
it appears that superexchange in these systems is suffi- 
ciently strong to destabilise the homogeneous ferromag- 
netic state at i < 1 [38]. The abrupt wall picture as 
discussed here is then inapplicable (see Sect. IV be- 
low). Wc note, however, that this might not be the 
case for other lightly doped magnetic semiconductors or 
semimetals. Ferromagnetic semiconductors such as Eu- 
doped EuS and EuO have relatively high values of Curie 
temperature Tc (Ref. [39]), presumably originating from 
a strong ferromagnetic superexchange, J < 0. In this 
case, even in a lightly-doped sample ferromagnetism is 
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due mostly to superexchange (rather than to double ex- 
change) and one expects that the domain walls will be of 
Bloch type, like in conventional Heisenberg ferromagnets. 
However, other magnetic semiconductors such as EuSe 
become ferromagnetic only upon small electron doping 
[40] . In this case of small positive J, domain walls may 
in fact be abrupt. This also may be the case in a ferro- 
magnetic semimetal EuB 6 (Ref. [41]). It would therefore 
be most interesting to study experimentally the domain 
wall structure (in particular, the effect of domain walls 
on the transport properties) in the ferromagnetic films of 
these compounds. 

Throughout our calculation, we neglected the effects of 
chemical disorder which can lead to localisation of elec- 
tron states. We note that the overall profile of carrier 
wave functions does not directly affect the properties of 
an abrupt wall. The assumption essential for our ap- 
proach is that the electron wave function can be locally 
approximated by an energy eigenfunction of the clean 
case [42] with the same energy. This is valid provided 
that the localisation length is much larger than the in- 
verse Fermi momentum; the latter condition is expected 
to be satisfied in manganates within the metallic regime, 
as well as in the doped magnetic semiconductors and 
scmimetals discussed above. 



IV. PHASE SEPARATION AND STRIPE WALLS 

Phase separation is a phenomenon which commonly 
occurs in the CMR manganese oxides [20,21]. Although 
direct evidence is lacking, it appears likely that the films 
studied in Rcfs. [5-7] are in fact phase-separated. It is 
therefore important to consider the effect of phase sepa- 
ration on magnetic domain wall structure in double ex- 
change ferromagnets. 

Let us first suppose that the values of parameters of 
the system (that is, carrier density x, superexchange J, 
Hund's rule coupling strength Jh) lie within the stability 
region of the uniform ferromagnetic phase. The thermo- 
dynamic potential is then given by £Ifm = E + 2 J — fix, 
with the value of fi = Hfm{x) determined by the uni- 
form conduction electron density x. The electron charge 
density, -ex, is compensated by the combined charge 
of magnetic and non-magnetic ions, resulting in electric 
neutrality of the system on the microscopic level [43] . As 
the values of parameters are varied {e.g., either the value 
of J is increased or that of Jh is decreased), the sys- 
tem eventually becomes unstable with respect to phase 
separation into ferromagnetic phase and another phase 
which we will call antifcrromagnetic [44]. In the absence 
of Coulomb interaction, this occurs when the thermody- 
namic potentials of the two phases become equal to each 
other: £Ifm{hfm{x)) = VIafm{ij-fm{x)). At this point, 
it becomes energetically advantageous to create islands of 



the antifcrromagnetic phase within the bulk of ferromag- 
nct. Since there is a finite energy cost W associated with 
a unit length of the boundary between the two phases, 
such an island should contain a large number of sites in 
order to reduce the boundary energy per antiferromag- 
nctic site; as long as this is the case, the area occupied 
by the antifcrromagnetic phase can be arbitrarily small 
relative to the total size of the system, so that the carrier 
density x within the ferromagnetic area and hence the 
value of chemical potential hfm{x) remain unchanged. 

Structure of the boundaries between different phases 
has been studied by the present writer in Ref. [22]. It 
was found that at least in some cases these boundaries are 
abrupt; it appears plausible that this property is rather 
generic. We note that the energy and charge of an abrupt 
interphase boundary can be evaluated using the approach 
applied in Sect. Ill above to the study of abrupt domain 
wall. A boundary between ferro- and antiferromagnetic 
areas can be perfectly abrupt only if it runs parallel to 
certain lattice directions [22]. It is therefore likely that 
within a large region of parameter values, the emerging 
islands of antifcrromagnetic phase will have a square (or 
diamond) shape. Apart from one case discussed towards 
the end of this section, the latter feature is unimportant 
for the rather qualitative discussion below. We will there- 
fore assume that the islands are circular, which would 
correspond to the boundary energy W independent on 
direction. 

While the chemical potential fi = [1fm{x) is constant 
across the sample, the carrier density within the island, 
xafm, is different from the nominal value x. We note 
that phase separation consists precisely in a redistribu- 
tion of the carriers with a simultaneous change in mag- 
netic ordering, and would not be possible had the re- 
quirement of constant carrier density been enforced on 
the microscopic level. The island is therefore electrically 
charged, and it is imperative to take into account the ef- 
fects of electrostatic Coulomb interaction and screening 
on phase separation. 

In a thin film, the inverse Dcbyc Huckel screening ra- 
dius is given by [45,46] (see also Appendix B): 



«= = — , e = + 6 d2) ■ 



(37) 



Here vq is the value of carrier density of states at the 
Fermi level and edi,e<i2 are dielectric constants of the me- 
dia on both sides of the conducting layer. In the 3D 
case, which is discussed in more detail in Appendix B, 
K (3D) = 47re 2 i' /e d , where e d is the dielectric constant 
of the double exchange magnet itself. If the size of the 
island was large in comparison with Debye-Hiickel ra- 
dius, R 3> k _1 , screening within the island would have 
restored the carrier density to its nominal value x (and 
charge density to zero) . In the case when there is no con- 
duction band in the bulk antifcrromagnetic phase {e.g., 
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when xafm(p) equals either or 1), the presence of elec- 
tric potential ip (which in this case is strongly position- 
dependent) would shift the carrier band within the is- 
land either upwards or downwards. This in turn will 
ultimately give rise to a Fermi surface, screening, and 
restoration of the carrier density to its nominal value on 
the length scale of However, as explained above, 

when the value of density is fixed no phase separation 
is possible. We therefore conclude that formation of 
an island can be energetically favourable only as long 
as [47] R ~ kT 1 . We will assume for simplicity that 
1 <C R <C that is, that the carrier density within 

the island is uniform and equal to the bulk value of 
xafm(p)- This obviously includes an assumption that 
Debye-Hiickel radius is large on the atomic length scale, 
k -C 1. The latter is not unphysical, in view of rela- 
tively large dielectric constants 6d reported for the man- 
ganates [48] and of suppression of carrier density of states 
at the Fermi level vq found in the x-ray absorption and 
angle-resolved photoemission measurements [49]. For a 
thin film, the situation also depends on the choice of the 
substrate, as discussed in more detail in the end of this 
section. 

With these assumptions, the change in thermodynamic 
potential Q associated with a creation of a single circular 
antifcrromagnetic island in a 2D system can be evaluated 
as [50] 



X t = nR 2 (n AFM - n FM ) + 2nRW + X - J d 2 rp{r)p{r) 



I d 2 r / de(e — p)v n . 

J{FM) J u. 



(38) 



Here, the first two terms represent the bulk and boundary 
contributions, the third term is the electrostatic Coulomb 
energy, and the last term is the kinetic energy cost of 
re-distributing electrons in the ferromagnet, caused by 
the shift of electrochemical potential p'(f) = p + e<p{r) 
(that is, the shift of band energies due to the presence of 
electric field within the screening cloud). Charge density 
p(r) equals pafm = —&{xafm — x) within the island, 
and — eSx(r) outside, where Sx is the change of electron 
density in the screening cloud. The last term in Eq. (38) 
can be re-written as 



e<p5xd 2 r 



ppd 2 r . 



I(FM) "0 z J{FM) z J(FM) 

This allows us to render Eq. (38) in the form 
X t = irR 2 (n AF M - SIfm) + 2ttRW + \ [ d 

2 J(AFM) 



rpip 



(39) 



where the integration in the last term is carried out over 
the area of the island. Evaluating the potential ip to lead- 
ing order in l/{nR) <C 1 [see Appendix B, Eqs. (Bll 
B12)], we obtain 



Xi = irR 2 (n AF M - n FM ) + 2ttRW + 



8irp 2 AFM R 3 



3e 



(40) 



Creation of an island becomes energetically favourable 
once the minimum value of this expression drops below 
zero. This yields the following threshold condition for 
the phase separation to occur: 



Q.FM — &AFM > A 



(41) 



Ao, the discriminant of the 
= vanishes; the minimum 
then reached at R = 



is 



[at Qfm — &afm - 
cubic equation Xi (R) 
value, Xi(Ro) = 

{ZeWyi 2 /(2\ P AFM\)]. 

Let us now consider a domain wall in a phase-separated 
film. We note that in this case the antifcrromagnetic 
and ferromagnetic tendencies in the system are approx- 
imately balanced against each other; this greatly re- 
duces both the spin stiffness [which in turn determines 
the Bloch wall energy via Eq. (2)] and the energy of 
abrupt domain walls, Sf d . This point is illustrated by 
Fig. 6, representing the chemical potential dependence 
of spin stiffness (solid line) and abrupt wall energies 
(dashed and dashed-dotted lines) for a Jh — > oo sys- 
tem with the value of J = J(p) adjusted in such a way 
[51] that Qfm = &afm- The appropriate antiferromag- 
netic phase near the endpoints p = ±2 is characterized 
by the usual Neel {n,Tr} (G-antiferromagnetic) spin or- 
dering, whereas in the vicinity of quarter-filling, p = 0, 
the A-antiferromagnetic phase with the ordering vector 
{tt, 0} proves more advantageous. The plethora of pos- 
sible phases arising in the intermediate case (see Ref. 
[22]) are not considered here, and no value is plotted for 
DS and unless the phase separation into the ferro- 
magnetic and either G- or A-antiferromagnetic phases is 
possible. Comparing Fig. 6 with the J = case, plotted 
in Figs. 1 (a) and 4 (a), we find a drastic reduction of 
both spin stiffness and domain wall energies at the in- 
termediate doping values. In addition, the energies of 
abrupt domain walls are now of the same order of mag- 
nitude as spin stiffness, in a marked difference with the 
single-phase case considered earlier. 

We will first discuss the effect of Coulomb forces in 
the case when the value of Qfm — Qafm is just above 
the threshold, Eq. (41), so that the islands of antiferro- 
magnetic phase arising within each ferromagnetic domain 
are well separated from each other, and Eqs. (38-40) are 
valid. As discussed in Sect. Ill, the abrupt domain wall 
shuts the carrier hopping in the perpendicular direction, 
acting as a partition in the gas of conduction electrons. 
In the absence of Coulomb forces, the energy cost of cre- 
ating a stripe of antifcrromagnetic phase adjacent to the 
wall is therefore equal to — (SIfm — ^afm)<1 (where d is 
the stripe width) per unit length of the stripe, and does 
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not include any additional boundary contribution. This 
statement (which is equivalent to saying that the abrupt 
wall energy is equal to 2W per unit length) is exact in 
the Jh —> oo limit (see Appendix A and Fig. 10). It 
is also clear that it provides a reasonable estimate for 
the case of large but finite Jh; the details of situation 
at finite Jh will be addressed elsewhere. Thus, when 
f2 fm — &afm > 0, the energy of an abrupt domain wall 
can be further lowered by inserting alongside it a stripe 
of antifcrromagnctic phase (sec Fig. 7). The width of 
the stripe is determined by a trade-off between the bulk 
and Coulomb energies, i.e., by minimising the energy of 
a stripe domain wall per unit length, 

d 2 

S s (d) = 2W - {Q FM - n AFM )d - —p 2 AFM \nKd (42) 

[see Appendix B, Eqs.(B6-B8)]. Since the antifcrro- 
magnetic stripe separates two ferromagnetic domains 
with antiparallel directions of magnetisation, the spins 
at the two edges of the stripe must point in the oppo- 
site directions. For the stripe of A-antifcrromagnet (G- 
antifcrromagnet) parallel to a lattice direction (lattice 
diagonal), this means that the number d (the number 
\/2d) must be odd [52]; similar conditions should hold 
for other phases. Since we assumed that the value of d is 
sufficiently large, d 3> 1, these requirements do not affect 
our estimates. Assuming that Q F m — &afm = AQoi we 
find 



d s = 



S s (d )^2W- 



3 InKcL 



(43) 



Eqs. (42-43) are valid to leading order in nd <C 1; 
even though hiKd s is thus large, the relatively large co- 
efficient of 16/3 in the second term of Eq. (43) allows 
for a significant reduction of domain wall energy due to 
the presence of a stripe of antiferromagnetic phase. It is 
not impossible that this reduction can make the quantity 
S s lower than the Bloch wall energy Sb, provided that 
the easy-axis anisotropy constant K is sufficiently large. 
The domain walls would then have a stripe structure, and 
would strongly interfere with the transport properties of 
the system. However, the exact values of quantities k 
and W in Eq. (43) arc not known, and it is not clear 
whether this situation can be realized experimentally. 

More importantly, Eq. (43) [and its 3D analogue, Eq. 
(B15)] refer to the case when phase separation is just 
beginning, with the islands of antiferromagnetic phase 
separated by large areas of a ferromagnet. Indeed, our 
derivation relied on an assumption that the screening 
clouds formed around different antiferromagnetic islands 
do not overlap, that is, that the inter-island distance is 
much larger than the screening radius. The size of each 
island, on the other hand, is much smaller than k _1 , 
so only a small part of the net sample area is occupied 



by the antiferromagnetic phase, making phase separa- 
tion difficult to detect. The available experimental data 
on phase separation in the CMR compounds [21], on the 
other hand, correspond to the case when a substantial 
part of the sample reverts to a non-ferromagnetic phase. 
Within the context of phase separation mechanism con- 
sidered here this is only possible when neither the size of 
antiferromagnetic islands (or stripes) nor the inter-island 
distance is larger than Debye-Hiickel radius. Below we 
will consider the case when screening is negligible (that 
is, when the inter-island distance is much smaller than 
k _1 ). Since k is expected to be small (see above), this 
is not unrealistic; moreover, the results are expected to 
provide a reasonable estimate for the case of intermediate 
screening strength as well. 

The ferro- and antiferromagnetic phases arc then char- 
acterized by uniform values of electron densities x F m and 
xafm and charge densities, p F u = —e(x F M — x) and 
Pafm — ~e(xAFM — x). The numbers of sites occupied 
by ferro- and antiferromagnetic phases, 



N F m — 



N 

TTs 



Nafm = 



NS 
T+6 



£ _ Pfm 



PAFM 



(44) 



(where N is the total number of sites in the system) 
are self-adjusted in such a way that the values of bulk 
thermodynamic potentials of the two phases, Ofm and 
Qafm, are close to each other. Therefore our observa- 
tion that both spin stiffness DS and abrupt wall energies 
are significantly reduced and are of the same order of 
magnitude (see above and Fig. 6) remains applicable. 

It is expected that the value of parameter S can be 
determined experimentally. 

We are interested in the situation when within each fer- 
romagnetic domain the poorly-conducting antiferromag- 
netic phase forms disconnected droplets (so that metallic 
conductance through the connected ferromagnetic area 
is still possible), and we will again assume that these 
droplets are circular in shape. The number of droplets 
in the sample is then Nafm/(kR 2 ) (where R is the ra- 
dius of a droplet), and thermodynamic potential of the 
phase-separated system is given by 



FM 



n 



AFM<- 



1 



5 



irR 2 l + 5 



(2ttRW + E 1 ) 
(45) 



per site, where E\ is the Coulomb energy of a single 
droplet. This term cannot be evaluated rigorously; in 
order to estimate it, we calculate the energy of Coulomb 
interaction within the so-called Wigncr cell, composed of 
the droplet and a surrounding ring R < r < R' (where 
r is the distance from the centre of the droplet) of the 
ferromagnetic phase. The value of R 1 = R[(l + 5)/5] 1 / 2 
is chosen in such a way that the combined charge of the 
droplet and the ring vanishes. It should be emphasised 
that unless 6 is small, 5 <C 1, this procedure, which has 
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been used to treat a similar problem earlier [53], is not 
exact [54]: even though the electrostatic potential of a 
Wigner cell falls off rapidly with distance, ip(r) oc r~ 3 , 
it does not vanish outside the cell. In addition, dif- 
ferent Wigner cells overlap with each other. Thus, by 
using this approach we essentially replace the Coulomb 
force with some model interaction, which however cap- 
tures the essential features of the original problem as long 
as the value of 8 is not too large (see below). We find 
E x = %Tx{B!fRp\ FM A x (b)bl{U) [see Appendix B, Eq. 
(B17)], where for small values of 8 <C 1 the function 
Ai(S) is equal to 1. The thermodynamic potential of the 
droplet phase, Eq. (45), has to be minimised with respect 
to the droplet radius R, yielding 



VLfm + QafmS , 8\pfm Ai(5)W 



1 + 8 



VTT6 



3e 



(46) 



Another possible geometry of phase separation is repre- 
sented by the stripe phase [shown in Fig. 8 (a)], formed 
by the parallel antiferromagnetic stripes of width d em- 
bedded into the ferromagnetic background. The thermo- 
dynamic potential of the stripe phase is given by 



n 2 (d) 



FM 



FM<~ 



1 8 



1 



dl + 8 



(2W + E 2 ), (47) 



Within the Wigner-cell approximation, Coulomb energy 
per unit length of a single stripe, E 2 , is calculated 
by subdividing the sample into the "Wigner stripes" 
of width d' = d(l + 8)/ '8 [see Fig. 8 (a)]. We find 
E 2 = -(ddf p\ FM A 2 (8)8/e)\n8 with A 2 {8 -> 0) = 1 [see 
Appendix B, Eq. (B19)]. Minimising the value of Q 2 
with respect to d, we obtain 



d = 



1 



2We 



n 2 - Oi 



\Pafm\ V {l + 8)A 2 (8)\\n8\ ' 

^ [2^2^2(5)111151 - 8 



(48) 




MS) 



The latter quantity is positive for all values of 8 be- 
tween and 1, indicating that within this model ap- 
proach, the droplet phase is always preferred (see below). 
The formation of a stripe domain wall in the droplet 
phase involves re-arranging spins within a Wigner stripe 
of width d' s = d s (l + 8)/8 into the stripe phase [see Fig. 8 
(b)], that is, forming a single stripe of antiferromagnetic 
phase [of width d s {8)\ flanked by two stripes of ferromag- 
nct. The net area occupied by antiferromagnetic phase, 
Nafm, is conserved, as is the overall electric neutral- 
ity. Minimising the stripe wall energy per unit length, 
S s = (Q 2 (d s ) — Qi)d' s , with respect to d s , we find 



WlA x 8 



A 2 {8)\p AFM ln8\]] 3(1 + 8) 
S s = 2WB{6) , B{8) = 1 - | 



MS) 



3A 2 (S)\\n8\ ' 



(49) 
(50) 



The ratio, B(8), of the stripe wall energy S s to the abrupt 
wall energy, 2W, is plotted in Fig. 9 (solid line). We 
see that the inclusion of an antiferromagnetic stripe can 
lower the energy of an abrupt wall by a factor of 4. Since 
the spin stiffness DS is of the same order as the abrupt 
wall energy (see Fig. 6), the stripe wall energy can be 
lower than the Bloch wall energy Sb already at a mod- 
erate value of anisotropy, K <~ DS/QA [cf. Eq. (2)]. 

Within the Wigner-cell approach for circular droplets 
the other droplet phase, with the ferromagnetic droplets 
in the antiferromagnetic background, becomes preferred 
at 8 > 1 (cf. Ref. [53]). While this transition might 
give rise to new possible domain wall structures near 8 = 
1, this is not expected to be physically relevant due to 
the intrinsic limitations of the method. As the value 
of 8 increases towards unity, the Wigner cell estimate 
for Coulomb energy becomes progressively less reliable 
due to decreasing separation between the droplets. It 
is perhaps even more important that the effects of the 
droplet shape can no longer be ignored. 

As mentioned above, it is likely that the optimal shape 
of antiferromagnetic droplets is square; this would be in 
line with earlier results for double exchange model [22], 
as well as with the numerical results for phase separation 
in other similar systems [55]. In order to calculate the 
energy of the square-droplet phase at small 8, one can 
still use the Wigner-cell approach. Due to the increase in 
the droplet boundary energy, the combined Coulomb and 
boundary contribution to the thermodynamic potential 
of the droplet phase [the last term in Eq. (46)] increases 
by some 6 %. This in turn leads to a noticeable decrease 
in the quantity B(8) (dotted line in Fig. 9). 

As the value of 8 increases, the Wigner-cell method 
becomes completely unsuitable for the analysis of the 
square-droplet phase. Indeed, at 8 = 1 (that is, at 
Nfm — Nafm) the square-droplet phase corresponds to 
a checkerboard arrangement of equal ferro- and antiferro- 
magnetic squares, which has nothing in common with the 
Wigner cell picture (cf. Fig. 8) . It is therefore clear that 
thermodynamic potential of the square-droplet phase at 
sufficiently large 8 is well above the value given by Eq. 
(46). Accordingly, Eq. (50) significantly over-estimates 
the value of B{8) and hence the stripe wall energy, S s . 
While leaving this subject for future investigation, we 
note that it is entirely possible that at a certain value 
of 8 = 8 C < 1 thermodynamic potential of the square- 
droplet phase exceeds that of the stripe phase, Eq. (47) . 
The quantities B(S) and S s will vanish at this point [56], 
as exemplified schematically by the dashed-dotted line in 
Fig. 9. In this case, for any finite value of the anisotropy 
constant K > and sufficiently small 8 C — 8 > 0, mag- 
netic domain walls within the conducting phase would 
have stripe (as opposed to Bloch-like) structure. 

We close with a brief comment on the applicability of 
our analysis to the finite-thickness films. The results of 
Sections II and III for the spin stiffness and abrupt wall 
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energies are valid only as long as the carrier velocity com- 
ponent perpendicular to the film is negligible. However, 
our conclusion that on the brink of phase separation DS 
and the abrupt wall energies are generally of the same 
order of magnitude (as illustrated by Fig. 6) is likely to 
remain valid in 3D as well. Our assumption that screen- 
ing has the two-dimensional character is valid as long as 
the film thickness remains small in comparison with the 
two-dimensional Debye-Huckel radius, The film is 

then thin from the viewpoint of electrostatics [cf. Eq. 
(B5)], that is, there is no electric field in the perpendic- 
ular direction within the film [57]. The latter holds pro- 
vided that the film itself is homogeneous in this direction, 
i.e., that characteristic length scale of a phase-separated 
sample (the droplet radius, R ~ (IW) 1 / 2 /\pafm\), is 
larger than the film thickness. Given the typical experi- 
mental observations [21] that phase separation occurs on 
the scale of at least 50-100 nm, this last condition is not 
particularly restrictive. 

The Debye radius can be roughly estimated by as- 
suming that uq is of the order of inverse bandwidth 
(At ~ 5eV) divided by the unit cell area (~ 0.15nm 2 ). 
Taking in Eq. (37) 6^2 = 1 (dielectric constant of the 
air), we then find kT 1 ~ (e<ji + 1) • 0.08nm. The sub- 
strate used in the measurements of Ref. [6], lanthanum 
aluminate, has the dielectric constant [58] of e<ji w 24, 
resulting in k _1 ~ 2nm. It is therefore tempting to as- 
sociate the reported domain wall resistance [6] (large for 
the thinnest Pr 2 /3Sr 1 / 3 Mn03 films of 4 nm, vanishing 
for films thicker than 20 nm), which is observable below 
the Curie temperature, Tc ~ 130-RT, with the stripe walls 
which arise only as long as the the thickness of conduct- 
ing layer (which is presumably somewhat thinner that the 
film itself) is not large [59] in comparison with k . We 
note that the film thickness required for the lattice pe- 
riods (and hence the anisotropy constant, K, and Bloch 
wall energy) to relax to their bulk values is of the order of 
500 nm (cf. Ref. [7]). Thus, our suggestion provides an 
(hitherto lacking) interpretation for the disappearance of 
domain wall resistivity in the films thicker than only 20 
nm. 

The experiments of Ref. [5], on the other hand, were 
performed with (ferroelectric) strontium titanate sub- 
strate, with [60] Cdi « 1200 at T = l\QK, which yields 
kT 1 <~ lOOnm. The Lao.7Cao.3Mn03 (with the Curie 
temperature Tc = 250K) film [5] was 200 nm thick, 
and the domain wall contribution was observable below 
T = 11QK. Given the strong dependence of e<ji on tem- 
perature (e<ji w 24, 000 at low T, e<<i w 300 at room tem- 
perature), it appears plausible that domain walls have 
stripe structure at low temperatures, when the film thick- 
ness is not large in comparison with Furthermore, 
it is not unlikely that the above-mentioned transition at 
T w 110-ftT is due to the violation of this condition at 
larger T, and associated change of the domain wall struc- 
ture [61]. We emphasize that this discussion is specula- 



tive at best, as we make no attempt to adequately de- 
scribe the crossover between two- and three-dimensional 
screening nor to take into account the peculiar geometry 
of the sample used in Ref. [5] . 

It appears that stripe wall formation is in principle 
also possible in the opposite limiting case of a bulk 3D 
material, although the Wigncr-cell estimates given in Ap- 
pendix B suggest that somewhat higher values of 8 are 
required. The values of Debye radius, K~^ D y and dielec- 
tric constant, e^, of doped manganates are, however, not 
known, and, crucially, very small values of anisotropy 
make the Bloch wall energy very low. It is therefore ex- 
pected that in the 3D case the energy of Bloch wall is 
generally lower than that of a stripe wall, in agreement 
with the fact that no observable domain wall contribution 
to resistivity was reported for the manganate crystals. 

V. DISCUSSION 

In this article we showed that there are at least three 
different possible types of structure of a ferromagnetic 
domain wall, all of which can be realized within the dou- 
ble exchange model. The energies and charges of Bloch, 
abrupt, and stripe domain walls are also different, as are 
their anticipated contributions to the resistance and mag- 
netoresistance of the sample. The conventional, weakly 
charged Bloch walls (Sect. II), which generally arise in 
single-phase samples, become unstable at low carrier den- 
sities, when the abrupt walls (Sect. Ill) are preferred. 
For a phase-separated system, however, there is a region 
of parameter values when the domain walls acquire stripe 
structure (Sect. IV), characterized by a stripe of antifcr- 
romagnetic phase separating the two domains. 

It is not yet known whether all three types of wall 
can occur in the CMR manganate compounds. As fol- 
lows from the discussion in Sect. Ill, abrupt walls are 
expected to arise at low values of electron doping [62], 
x <C 1, provided that the homogeneous ferromagnetic 
phase remains thermodynamically stable. We are not 
aware of any measurements of the domain wall contri- 
bution to transport in this regime, and it is not clear 
whether such a situation (which also requires the value 
of direct supcrcxchangc J to be extremely small) can be 
realized in the manganates (however, see the end of Sect. 
Ill for a discussion of other compounds). As for the in- 
termediate doping values, it appears that domain walls 
can have either stripe or Bloch structure. 

The effect of Bloch walls on the charge transport prop- 
erties of a double exchange ferromagnet has been dis- 
cussed theoretically [63]. The results are consistent with 
simpler estimates [5-7] suggesting that for a realistic 
value of Ib and at an intermediate doping level, carrier 
scattering off the Bloch wall cannot possibly account for 
a measured domain wall contribution to the resistivity 
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of the systems studied in Refs. [5-7]. Measurable do- 
main wall contributions to the transport properties of 
the CMR manganates, reported in other studies known 
to us, are attributable to the grain boundary effects in 
polycrystalline films [9,11,64]. In this case, the magnetic 
structure [10,11] of a domain wall arising at a substrate 
grain boundary is largely determined by underlying lat- 
tice defects [65]. It is anticipated that this also holds for 
the magnetic pattern appearing in a strained film at the 
boundary of a heavy- ion irradiated region [66] . We note 
that the effects of lattice irregularities of any type are not 
included in the present theoretical treatment. 

Our results suggest that magnetic domain walls in 
monocrystals or epitaxial films of CMR manganates at 
the intermediate doping levels generally have Bloch-like 
structure, with a notable exception of certain strained 
films similar to those used in Refs. [5-7]. Regarding the 
latter case, we expect that domain walls may in fact be 
the stripe walls introduced in Sect IV above. This sugges- 
tion is corroborated by especially strong effect reported in 
Ref . [6] , which shows that domain walls give a dominant 
contribution to the resistivity of a thin Pr 2 /3Sr 1 / 3 Mn03 
film at low temperatures. The connexion between do- 
main wall resistivity and dielectric properties of the sub- 
strate, discussed in the end of Sect. IV, appears to 
lend further support to the stripe wall scenario. The 
stripe walls appear likely to arise in this case due to the 
strain-induced increase of easy-axis anisotropy constant 
K (which in turn increases the Bloch wall energy), and 
also to phase separation which makes formation of the 
stripe walls possible. While it is not clear whether phase 
separation does occur in the samples used in Refs. [5-7], 
this would be rather plausible given that phase separa- 
tion is commonly observed in both manganate crystals 
and films [21]. We suggest that further measurements (e. 
g., scanning tunnelling spectroscopy) need to be carried 
out to clarify whether these samples are in fact phase- 
separated. On the other hand, domain wall properties 
(including possible domain wall contribution to the re- 
sistivity) of those CMR films which are known to phase- 
separate [21,67] should also be investigated. Synthesis of 
electron-doped manganate films, if technologically pos- 
sible, may represent a promising new direction [38]. We 
note that magnetic domain walls appear only when a sub- 
stantial fraction of the film is in the ferromagnetic state, 
allowing for a low-field metallic conduction. 

In the present article, we did not quantitatively ad- 
dress the problem of conduction across a domain wall of 
either type. The available theoretical estimates of do- 
main wall conductance (Ref. [63] for Bloch walls, Ref. 
[16] for abrupt wall) are incomplete in that the Coulomb 
interaction between the carrier and the (charged) domain 
wall is not taken into account. As for the stripe walls, 
the issue of magnetotransport in this case has yet to be 
treated theoretically, although it is clear that stripe wall 
contribution to resistivity is much larger than that of 



either Bloch or abrupt walls. In the presence of stripe 
domain walls, magnetoresistance will be affected by the 
change of their structure under a magnetic field, which 
is likely to include a field-driven transition from stripe 
to Bloch walls. It is therefore expected that the depen- 
dence of the domain wall contribution to resistivity on 
the magnitude of applied in-plane field can be different 
for the Bloch and stripe cases (smooth decrease for Bloch 
walls, as opposed to possibly step-like features for stripe 
walls, as seen in Ref. [5]). 

Magnetotransport studies are not the only way to in- 
vestigate the properties of magnetic domain walls. Direct 
probes of charge and spin structure of domain walls are 
possible in principle (cf. Ref. [68]), but have not yet been 
performed for the manganates. However, Fresnel imag- 
ing measurements on a thin Lao.7Cao.3Mn03 film were 
reported recently [8]. Domain walls were found to re- 
tain a finite width of the order of 40 nm, in apparent 
agreement with Eq. (2) for Bloch walls. We hope that 
domain wall widths in the strained films studied in Refs. 
[5-7] will also be measured in the near future. It would 
be most interesting to try to relate these to the band 
structure, magnetic, and electrostatic properties of the 
corresponding compounds and to check the agreement 
with the estimates (43) and (49) for the stripe walls. 
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APPENDIX A: DERIVATION AND ANALYSIS 
OF EQS. (33-34). 

The key step in the calculation of the spectral shift 
function, Eq. (32), is the diagonalisation of perturbation 
operator, Vk x [see Eqs. (24) and (25)]. Its eigenvalues Ai 
and the corresponding fermionic operators dj are given 
by 

A 1 = A 2 = -A 3 = -A 4 = -^=^ 2 , (AI) 

A 5 = A 6 = -A 7 = -A 8 = — ^(1-V), (A2) 

V2 

and 

1 /2 

8 (2 ± V2) 01,3 = (4 - i/> 2 )(d-i T + d 2] ) + [(4 ± 2V2)V> + 
+ (4 ± V2> 2 ](doT + dif) + (1 ± V2)(4 - i> 2 ){d- n - d 2l ) + 
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+ [=F2v / 2'0 — (2 ± 3\/2)?/> 2 ](do|. — dii) , (A3)diagonalise the perturbation [see Eq. (25)] in the case of 

/ xi/2 a vertical wall as well. 

8 (2 T V2J a 2 , 4 = -(4-?/> 2 )(cLi T -d 2T ) + [(4 T 2V2)V> + We note that in the expressions for both H kl and V kl 

+ (4T^ ( * r -* t)+ (-i ± v5)(4-#»)( i . lt+ A t ) + - jrcaft^jEK is 

(A4)term on the r. h. s. of Eq. (A9), regardless of whether 



+ [±2\/2> - (2 T 3V2)^ 2 ](d oi + d n ) . 



/ r-\ V 2 r . , ,,,,, , . , ^- the domain wall is present, is to shift all of the energy 

8 2=fV2 a 5 7 =-2V2 (ip+ip 2 ) (d-it+ d 2 T + V2=f1)x , , , ' ~ 

V J ' 1 l; 1 T ' levels by - cos fci . Thus, Eqs. (26-27) are now replaced 



:[4± 2(x/2±l)V- -^ 2 ](rfo T +rfiT) ± ( 4 T 2>/2)(V> + V 2 ) x 



by 



x(d_ u -d 2i ) + [±4-2(V2 T l)V'T 2^ 2 ](do| -d u ), (A5) 

1 /2 

8(2± V2) a 6 ,8 = 2V2(^ + V7 2 )(d_ 1T -d 2T ) + (V2±l)x 

x[4T2(V2Tl)^-^ 2 ](d 0T -rfiT)±(4±2V2)(^ + ^ 2 ) x 

x (d_ u + dai) + [T4 - 2(V2 ± 1)V ± ^ 2 ](d | + d u ) . (A6) 



= 



<50 



£1,2 



J ~2~tt~ J de 2<We 2 V(e 2 + d) = 
= |f^|rfe 2 |(e 2 )/(e 2 + e 1 ), (A10) 
= -cosfci, 2 , vtoti{c-i) = L 2 /(tt\J 1 — e|) . 



These expressions are then used to form the matrix ele- 
ments My (see Eq. (29)), for example 

M n = ^V T 2 - 1) - ^(4^ - 1) + (2 - V2) x 
x(l + 3£ T - 4£ T 3 )/ T + (2 + V2)(l - 3E L + AE^)I L + Oty), 



Here, L\ and L 2 are the dimensions of the sample, and 
£(e) is the spectral shift function of the corresponding 
ID problem. It is evaluated as 7r£(e) = — ArgDct(<5ij — 
MijAj), with 



/„ = 



1 



dkn 



8ttQ J E a + cos k y 



iO 



(A7) 



Mii = £ / 



a=T,i" 



rffc 2 <0|a«|fcg><fcg|at|Q) 
27T £7 a + cos fc 2 — iO 



(AH) 



[see Eq. (30)]. Note that, owing to the symmetry proper- 
ties of the operators a,, the quantities My vanish unless 
both indexes i and j are either odd or even. Hence the 
8x8 determinant on the r. h. s. of Eq. (28) reduces to 
a product of two 4x4 determinants. After some algebra, 
one obtains expression (32), which has to be substituted 
into Eqs. (33) and (35). 

In the case of a vertical wall we choose the co-ordinate 
axes r\ and r 2 along the lattice directions with the r 2 
axis perpendicular to the wall. After the Fourier trans- 
formation, 

d Q (ri,r 2 )=7V- 1 / 4 ^e ikiri d a (/c 1 ,r 2 ), (A8) 
fei 

we find that the unperturbed Hamiltonian has the form 
n = J2\' H k 1 -cosk 1 J2dl(k 1 ,r 2 )d a (k 1 ,r 2 ) \ , (A9) 

fei V r 2 ) 

and domain wall again results in a local perturbation, 
Hhi — * "Hki + V/jj . This perturbation is still illustrated 
by Fig. 3, although the intersite distance is now equal 
to unity, rather than to 1/V2. The operators Hk ± and 
Vfej have the same form as Hk x and Vk x [see Eqs. (23- 
24)], with the substitutions Q — > 1, d a (k x ,y + 1/V2) — * 
d a {k\,y + 1), d a (k x ,i/\/2) — » d a (ki,i). Hence the Eqs. 
(A1-A6) with the value of Q set to unity can be used to 



and Ef = E± + Jh = e. Taking also into account that 
the states |fc 2 ) are defined in a conventional way, |fc 2 ) = 
Y, T2 exp(-ifc 2 r 2 )<4(fci,r 2 )|0) [cf. Eq. (31)], we conclude 
that the value of My coincides with that of My , Eq. 
(29), calculated at Q = 1. Thus, £(e) = £(e,Q = 1) [see 
Eq. (32)], and Eqs. (34), (36) follow. 

The somewhat cumbersome expressions (32-36) be- 
come much simpler in the case of large Hund's rule cou- 
pling, J H — > 00. We then find [69] 



S d = 



1 iV]EZ h m 

= < — arems 



arccos- 

7T 2 



V2 1 

+£ - M [a; - 0(h)] 1 + 2^2^ + X + 



4J+ 



+ 



3tt 2 



14 



1 + 



Y 2 



S v — 



2tt 



2tt 



-arccos(|yu| 



-H[x - 9(fj,)} + 4 \ J + K - — 



1) 
3 



m 2 + |m 4 
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T 

e 



Yn 



V> 2 , 



(A12) 
2J + E- 

(A13) 



^^arccos^sgn^- -^(2; -(%)), (A14) 

(Tt, = — — arccos(|^| — l)sgn^i — e(x — 6(fi)) , (A15) 
2tt 
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where Yi and E are given by Eqs. (5-6). When de- 
riving the -0 = values of Sd and S v above, it is con- 
venient to use a calculation scheme somewhat different 
from that used in the finite- Jh case. Namely, the lo- 
cal perturbation we consider now (see Fig. 10) cor- 
responds to inverting the spins along a ID chain, not 
only shutting the carrier hopping but also introducing a 
single chain of an antifcrromagnctic phase. The latter 
circumstance can easily be accounted for by subtract- 
ing the difference of thermodynamic potentials between 
the antiferro- and ferromagnetic phases; the advantage 
of this method lies in a very simple form of spectral shift 
function, £(e,Q) = (l/2)sgne, corresponding to the per- 
turbation shown in Fig. 10. 

Another potentially important case when the integra- 
tion in Eqs. (33-36) can be carried out analytically is 
that of small electron densities, For any value of 

Jh 3> x, we obtain: 



APPENDIX B: STRIPE WALLS AND 
SCREENING 

In this Appendix, we are concerned predominantly 
with investigation of screening potentials and Coulomb 
energies of phase-separated states in a two-dimensional 
conductor. Let the values of dielectric constants of the 
media on both sides of conducting plane be e^i and td2- 
The method of images enables one to evaluate the po- 
tential of a point charge q located at a distance z from 
a plane separating the two dielectric media [14]. In the 
limit z — ► 0, we find that this potential at any point in 
space is given by qj (es), where s is the distance from the 
charge and e = (e d i + e d2 )/2. We therefore conclude that 
the electrostatic properties of this system are described 
by a Poisson's equation of the form 



eV 2 </? = -4np(r)6(z) . 



(Bl) 



+V2 |tt J h (l - yjl + 3^ x 2 + 4J + 2K J 4> 2 , (A16) 
3 V Jh 



+4^jTj H \ [ l-^l + —Jx 2 + J + Kjip 2 . (A17) 

It is instructive to note that expansion of these expres- 
sions in the case of Jh S> 1 shows that the leading 1/ Jh 
correction amounts to a renormalisation of the superex- 
change constant, J — > J + (ttx 2 )/(2Jh)- This is another 
illustration of an effective antifcrromagnetism being in- 
duced by a finite Hund's rule coupling, as discussed in 
Sect. II. We also see that at J = and to leading order 
in x <C 1, abrupt wall energy does not depend on orien- 
tation of the wall, S d ^ = si°^ , which is due to the carrier 
dispersion law being isotropic at low densities. 

Electric charges of unperturbed abrupt domain walls 
at x <C 1 , Jh are given by 



lx Jh+4 



ex^, (MS) 



a v w e\l ex 

7T 



i / Jh + 2 : 3/2 V7T 



J H + "-" J lI' (A19 » 



Finally, we also quote a 3D result for a vertical abrupt 
domain wall energy (per unit area) at x <C 1 and J H — > 

oo: 

9l/3o4/3 

S (3C) Rj ^_£_ 7r 5/3 a ,4/3_ 2J) ^ = (K (A2Q) 



Here, r = {x, y} is the 2D radius-vector in the plane, 
z axis is perpendicular to the conductor, and V is the 
usual 3D gradient. It is therefore only the effective di- 
electric constant, e, that will affect the values of physical 
quantities in this case [cf. Eq. (37)]. 

We begin with evaluating the potential of a charged 
string within the film. Assuming that the string coincides 
with the x axis, we re-write Eq. (Bl) as 

eVV(y, z) - [47reV(», *H - A™8{y)]5{z) . (B2) 

Here, a is the linear charge density of the string, and the 
first term on the r. h. s. accounts for a screening charge 
arising from the band energy shift by the electrostatic en- 
ergy, — &p. This is a standard Thomas-Fermi treatment 
of screening, valid in the long-wavelength limit. Upon 
one-dimensional Fourier transformation we obtain 



dz 2 



ki ) <p(k y , z) = ( 2KLp(k y , z) - 1^ j S(z). 



V 



(B3) 



Using the Green's function for Eq. (B3) (cf. Ref. [46]), 



g{k y ,z) 
we obtain 



00 dk z e ik * z 1 



-\kyZ\ 



(B4) 



,2irk 2 + k 2 2\k y 
<p(ky, z) = |F^ e " |fcyZ| {^f - 0)) • (B5) 



Hence at z = 0, y(k y ) = 2-Ka/[e(\k y \ + «)], and 

-2(7 



(cos Ky cmy + sin ny siny) 



(B6) 



where si and ci are sine and cosine integrals. At ny 3> 1, 
Eq. (B6) yields ip(y) « 2a/(eK 2 y 2 ). Along with the 1/r 3 
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decay of a screened point charge potential [45,46], this is 
in contrast with the well-known exponential behaviours 
found in 3D. 

Let us now consider the potential of an antiferromag- 
netic stripe of width d <C centred around the x axis. 
At y » d, it is given by Eq. (B6) with a = pAFAid, 
whereas at y -C kT 1 it should coincide with the un- 
screened potential of the stripe, 



%Pafm 



v{y) = - 
+(o +y) ln k + v\- d 



-- y) ln|--y|+ 



consi , 



(B7) 



At y 3> d, the latter expression takes the familiar form, 



2cr , 

= — ~ m |y| + const . 



(B8) 



The two regions, y 3> d and y <C k _1 , overlap, en- 
abling us to find the value of const in Eqs. (B7-B8), 
-2p AFM d(C + Ln/c)/e where C w 0.577 is the Euler's 
constant. Substituting Eq.(B7) into Eq. (39), we find 
the leading order expression for the Coulomb energy of 
the stripe [the last term in Eq. (42)]. It is also easy to 
estimate the Coulomb energy of a Bloch wall, 

1 f f°° 1 
--II Sx(y)\n\K(y - y')\Sx{y')dydy' ^ -a%\n\Kl B \ 

^ J J — CO ^ 

(B9) 

[see Eq. (18)], assuming that kIb <C 1. 

We note that screening affects the value of potential 
at 2 1 y | < d even in the limit of d <C kT 1 because the un- 
screened potential, Eq. (B8), diverges at y — > oo. This 
is also the case for the potential of a single antiferromag- 
netic layer of thickness d in a phase-separated sample in 
three dimensions, in which case we find, per unit area, 

= (BIO) 



L Jafm 



pipd 3 r 



j2 2 

7r " Pafm 
K (3D) e d 



K (3D) 



In the case of a single antifcrromagnctic disk in 2D, or an 
antifcrromagnetic sphere (ball) in 3D, the unscreened po- 
tential vanishes at large distances and the leading-order 
(in kR) term in the Coulomb energy does not depend 
on the screening radius. Indeed, for the 2D case at suf- 
ficiently small distances r <C k , the exact screened 
potential of a point charge irR 2 pafm (found in Refs. 
[45,46]) is to leading order given by nR 2 pafm / '(^r). For 
r>fi, this clearly matches the unscreened potential of 
a charged disk. Therefore screening does not affect the 
value of ip within the disk, which enters Eq. (39). Using 
the Green's function procedure similar to Eqs. (B3-B4) 
above, we obtain for an unscreened disk of radius R 



<p(r)= J ^)e ifcV 0, <p(k) = 



A-k 2 pafmR 
Ik 2 



Ji{kR) 
(Bll) 



(where r is the distance from the island centre, and J\ is 
Bcssel function), and 

(B12) 

leading to Eq. (40), whereas for a 3D sphere of radius R 
we readily find 



1 6tt 2 P 2 afm R 5 
afm 15ed 



pipd r 



(B13) 



Eq. (41), derived in Sect. Ill, holds for a thin film. With 
the help of Eqs. (BIO) and (B13), it is easy to obtain a 
similar phase-separation threshold condition for the 3D 
(bulk crystal) case, 



Qfm — Qafm > A 



1/3 



5e d 



(B14) 



Here, W( 3 d) is the energy per unit area of a 
fcrromagnetic-antifcrromagnetic boundary in three di- 
mensions, which can be approximated by half the value 
of the 3D abrupt wall energy (cf. Eq.(A20)). The en- 
ergy of stripe (layer) domain wall at the phase-separation 
threshold in 3D is then given by [cf. Eq. (43)] 



S (3D) = 2W 



(3D) 



per unit area, where 
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^Pafm 



1/3 



is the radius of antifcrromagnctic bubbles appearing im- 
mediately above the threshold, Eq. (B14). 

We now turn to the other regime of phase separation 
considered in Sect. IV. In this case, screening is negligi- 
ble and our estimates of Coulomb contributions to ther- 
modynamic potential are based on evaluating the elec- 
trostatic energy of a single unscreened Wigner cell. In 
the case of circular antiferromagnetic islands ("droplet 
phase" ) , the Fourier component of electric potential of a 
Wigner cell is given by [cf. Eq. (Bll)] 



Air l R' 
ek 2 R 



Pfm [RJi {kR 1 ) - R'Ji (kR)} ■ (B16) 



Momentum integration [cf. Eq. (B12)] then yields the 
expression for the Wigner cell energy E\, given in the 
text above Eq. (46), with 
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1 

2' 



2' ' 1 + 



(B17) 



16 



where 2 F\ is the hypergeometric function. 

For the stripe phase, we find the electric potential <p(y) 
of a single "Wigner stripe" , 



Pafm 



= 2yS\n 



d! + 2y 



+d'6ln 



d'-2y 
{d'f - 4y 2 



- 2y(l + <5)ln 



d + 2y 



d-2y 



c? 2 - 4y 2 



+ 



(B18) 



where y is the distance from the stripe centre. We note 
that at y ^> d' , the potential ip(y) decays as 1/y 2 . Eval- 
uating the electrostatic energy per unit length, E 2 — 

f-d ,2 /2 ( P(y)P(y}dy/2'i we obtain the expression given in 
the text following Eq. (47), where 



A 2 (8)ln5 = (l + 5)ln 



4(5(1 +8) 1 + 25 + 25' 
1 + 26 



25 



ln(l + 25) . 
(B19) 



It is also possible to evaluate the Coulomb energy of the 
stripe phase exactly, taking into account the interaction 
between different "Wigner stripes" . Numerical calcula- 
tion shows that this leads to an increase of the quantity 
A 2 by about 8% at 5 -> 1, and by only 3% at 5 = 0.17 
[the latter corresponds to the minimum of B{5) in Fig. 
9] , attesting to the relatively high accuracy of the Wigner 
cell method for the stripe phase even in 2D. 

Coulomb energies of droplet and layered phase- 
separated states in 3D were evaluated in Ref. [53]. The 
sum of boundary and Coulomb contributions to the ther- 
modynamic potential equals 



n(3D) 



n 



(R) = 



35 W, 



(3D) 



1 + 5 R 



An 
5e d 



R*p 2 AFM Af D \5)5, (B20) 
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for the droplet phase, and 



nf D \d) = 



25 W, 



(3D) 



1 + 5 d 



j2 2 

+ -^-d p AFM 



(B21) 
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for the layered phase. Minimising expression (B20) with 
respect to the radius R of spherical antiferromagnetic 
droplets, and then the 3D stripe wall energy per unit 
area, si 3D) = 2 D) (d s ) - Cl {3D) )d' s [where d' s = d s (l + 
6)/ 5], with respect to the antiferromagnetic layer thick- 
ness <i s , we find 



3 5 / 6 -y / 2^(4 3g) ) 1 / 6 ( W (m e d 
V/W^Pafm?' 3 \ 1 + 5 



1/3 



Sf D ^ = 2W im B im (5), B {3D) (8) = 1-3 A /-4 3D) «5. 

(B23) 



The ratio B( 3 d)(5) of the energies of stripe and abrupt 
walls in 3D is plotted in Fig. 9 (dashed line). We see 
that the stripe wall energy vanishes already within the 
Wigner-cell method as the value of 5 approaches 
0.47. 
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FIG. 1. Spin stiffness DS (a) and the coefficient C [see Eq. 
(19)] (b) vs. electron density x for J = and Jh — » oo (solid 
line), Jh = 8 (dashed line), and Jh = 4 (dashed-dotted line). 
Dotted lines correspond to a regime where the spin stiffness is 
still positive, D > 0, but the ferromagnetic phase is unstable 
with respect to phase separation. 



FIG. 2. Diagonal (a) and vertical (b) abrupt domain walls 
(dashed lines). 

FIG. 3. Schematic representation of a one-dimensional 
problem which arises in diagonal domain wall calculations, 
Eqs. (23-24). The intersite distance is equal to l/\/2, and 
the numbers are the same as the subscripts of the fermion 
operators in Eq. (24). Dashed arrows correspond to the per- 
turbed case, ip 0. 

FIG. 4. (a). Abrupt wall energies vs. x at J = 0. Solid 
(dashed) lines, top to bottom: diagonal wall energy, S^ (ver- 
tical wall energy, Si°^) for Jh — > oo, Jh = 8, and Jsr=4. For 
finite values of Jh, the lines end at the values of x correspond- 
ing to the sign change of spin stiffness, D. Immediately below 
these values, the ferromagnetic state is unstable with respect 
to phase separation (see Fig. 1). For K(x) = D(x)S/25, 
(i.e., Ib = 5) the quantities Z d>v are negative everywhere 
outside the low-doping regions x <§C 1 and 1 — x <§C 1, ex- 
cept for the case of Jh = 8, when Z d becomes positive for 
x > 0.83 (dotted line), (b). Abrupt wall charges in units 
of electron charge, |e|. Solid and dashed (dashed-dotted and 
dotted) lines represent a d and a v for Jh — » oo (Jh = 4). 
(c). Bloch wall energy Sb (solid line), abrupt vertical wall 
energy Si°^ (dashed line), and the quantity Z v (dotted line) 
vs. superexchange J. Anisotropy constant varies according 
to K(J) — DS/25. Conduction electron density and Hund's 
rule coupling strength are given by x — 0.55 and Jh = 4, 
respectively, and the system is unstable with respect to phase 
separation. 

FIG. 5. Bloch wall energy Sb (solid line), abrupt ver- 
tical wall energy S^ ' (dashed line), and the quantity Z v 
(dashed-dotted line) vs. electron density, x, in the low-density 
limit without superexchange (J=0). Hund's rule coupling 
is fixed at Jh = 0.1, while the anisotropy varies according 
to K(x) = DS/25, leading to a constant Bloch wall width, 
Ib = 5. The Bloch wall, however, becomes unstable at lower 
x (dotted line). 
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FIG. 6. Chemical potential dependence of spin stiffness DS 
(solid line) and the diagonal and vertical abrupt domain wall 
energies (dashed-dotted and dashed lines) for a Jh — * oo sys- 
tem on the brink of phase separation. The value of J is ad- 
justed in such a way that Q,fm = Qafm for any value of 
carrier density x. The nature of corresponding antiferromag- 
netic phases is discussed in the text. 

FIG. 7. Schematic representation of a stripe domain wall 
in a phase-separated double exchange magnet. The two fer- 
romagnetic domains with antiparallel directions of magneti- 
sation (arrows) are separated by a stripe of antiferromagnetic 
phase (shaded). In addition, unconnected islands of antifer- 
romagnetic phase are formed within each domain. 

FIG. 8. Stripe phase (a) and stripe domain wall within 
the droplet phase (b). The system is phase-separated into 
ferromagnetic (unshaded) and antiferromagnetic (shaded) re- 
gions with 8 w 0.4. The Wigner cell boundaries of stripe 
and droplet phases are shown with dashed and dotted lines 
respectively. The two connected ferromagnetic domains ex- 
tending to the left and to the right of the stripe wall in (b) 
are magnetised in the opposite directions (not shown). The 
width of antiferromagnetic stripes in (a) and (b) is given by 
Eqs. (48) and (49), respectively. 

FIG. 9. The ratio B of the energy of a stripe wall to that of 
an abrupt wall [see Eq. (50)] vs. the ratio S of antiferro- and 
ferromagnetic areas of the sample: solid line, droplet phase; 
dotted line, square-droplet phase at low S; dashed-dotted line, 
possible behaviour for the square-droplet phase at larger 5; 
dashed line, 3D result of Eq. (B23). 

FIG. 10. Local perturbation used in the calculation of 
abrupt domain wall energies at Jh — > oo (for the case of 
a diagonal wall). 
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